Stiff spinning torus of electromagnetic two-body motion 
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We study an orbit of the electromagnetic two-body problem that involves a fast (stiff) spinning 
motion about a circular orbit. We give a multiscale method of solution that solves for the fast 
timescale first. The solvability condition of the asymptotic expansion demands resonances for the 
fast dynamics. The stiff resonant tori are found precisely in the atomic magnitude and agree 
with many features of the Bohr atom in quantitative and qualitative detail; We calculate every first 
emission line of the first f 3 observable spectroscopic series of hydrogen within a few percent deviation. 
The resonant orbits have angular momenta that are approximate multiples of Planck's constant and 
the emitted frequencies are given by a difference of two linear eigenvalues. This Lorentz-invariant 
two-body dynamics exhibts the phenomenon of resonant dissipation, i.e., the metastable dynamics 
radiates the center-of-mass energy while the particles perform fast spinning oscillations of small 
amplitude about a circular orbit, a collective radiative recoil. 



We study a two-body orbit of classical electrodynam- 
ics involving a fast (stiff) spinning toroidal motion about 
a circular orbit. The nontrivial stiff dynamics is part of 
the full two-body motion with delay and it disappears 
in the Coulombian limit. Our equations of motion are 
prescribed by Dirac's electrodynamics of point charges 
with retarded-only fields , a fundamental physical the- 
ory. The existence of the fast delay-dynamics demands 
a multiscale asymptotic solution that solves for the fast 
timescale first. The multiscale method requires a res- 
onance between two fast frequencies because of Freed- 
holm's alternative Q , providing the natural coupling be- 
tween the fast and the slow timescales. The resonant 
orbits turn out to be precisely in the atomic magnitude; 
These resonant orbits are labelled by an integer and have 
angular momenta that are multiples of a basic angular 
momentum, that agrees well with Planck's constant. The 
emitted frequencies agree with the sharp lines of quantum 
electrodynamics (QED) for the circular lines of hydrogen 
with a few percent average deviation. The angular mo- 
mentum of the stiff spinning motion is of the order of the 
orbital angular momentum, again in qualitative agree- 
ment with QED. 

Dirac's 1938 work on the electrodynamics of point 
charges gave delay equations that were seldom studied. 
Among the few dynamics investigated, an early result 
of Eliezer 0, 0, @ revealed a surprising result (hence- 
forth called Eliezer 's theorem); An electron moving in 
the Coulomb field of an infinitely massive proton can 
never fall into the proton by radiating energy. The re- 
sult was generalized to tridimensional motions with self- 
interaction in a Coulomb field 0, Q , finding that only 
scattering states are possible. Since our model has the 
dynamics of Eliezer's theorem as the infinite-mass limit, 
a finite mass for the proton is essential for a physically 
meaningful dynamics; If the proton has a finite mass, 
there is no incrtial frame where it rests at all times, and 



this in turn causes delay because of the finite speed of 
light. A finite mass for the proton is what brings de- 
lay into the electromagnetic two-body dynamics, with 
its associated fast dynamics. We shall describe the two- 
body motion in terms of the familiar center- oj-mass co- 
ordinates and coordinates of relative separation, defined 
as the familiar coordinate-transformation that maps the 
two-body Kepler problem onto the one-body problem 
with a reduced mass plus a free-moving center of mass. 
We stress that in the present relativistic motion the 
Cartesian center-of-mass vector is not ignorable, and it 
represents three extra coupled degrees-of-freedom. We 
introduce the concept of resonant dissipation to general- 
ize the concept of non-ionizing orbits of Ref.0 . Reso- 
nant dissipation is the condition that both particles recoil 
together, i.e., the center-of-mass vector decelerates, while 
the particles perform the almost-circular orbit illustrated 
in Fig. 1, despite of the energy losses of the metastable 
center-of-mass dynamics. 

Here we use a circular orbit as the unperturbed solu- 
tion of our perturbation scheme. This circular orbit is 
a solution of the low- velocity Coulombian-limit obtained 
by ignoring the delay and its associated fast dynamics. 
The quantities for this uniform circular motion are given 
now |2j; The electron has mass m\ and travels a circle 
of radius r\ while the proton has mass mi and travels 
a circle of radius t-i- We henceforth use units where the 
electronic charge is e = — 1 and the speed of light is 
c = 1. The circular orbit is defined by the retardation 
angle 9 that one particle turns while the light emanating 
from the other particle reaches it (one light-cone distance 
away). The angular momentum l z of the circular orbit 
is given to first order by Zj = 9, and for atomic orbits 
9 is a small parameter of the order of the fine structure 
constant a — 1/137.036. For small values of 9 the fre- 
quency of angular rotation and the particle distance 
in light-cone rt, are given, to the leading order in 9, by 
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Kepler's law 



(1) 



where (i = m\mil{m\ + 777,2) is the reduced mass. The 
circular orbit is not an exact solution of the full two-body 
equations of Dirac's theoryQ] with a finite protonic mass 
and self-interaction; Substituting the circular orbit into 
the equations of motion of [l| , we find for the electronic 
motion an offending force opposite to the electronic ve- 
locity of magnitude 



(2) 



where we scaled the offending force by the magnitude 
of the attractive force between the particles, a scaling 
factor of \/t\ . Along the trajectory of the proton, the 
delayed Lienard-Wiechert attraction of the electron is the 
main offending force against the velocity a force of ex- 
actly the same magnitude of Eq. J2J) as calculated with 
the Page expansion of Ref.[8|. These offending forces 
cause dissipation along the unperturbed circular orbit. 
This radiative instability in the slow timescale is an in- 
complete picture of the dynamics though; it disregards 
the fast stiff motion that is present because of the delay 
0- Precisely because of the stiff delay dynamics, at very 
small deviations from circularity the nonlinear terms of 
the self-interaction force correct the leading term of Eq. 
(J2J , enough to balance the offending force of Eq. (J2J . In 
the following we explain in physical terms how the stiff 
torus of Fig.l is formed. 

We start by substituting the circular orbitplus a per- 
turbation into the equations of motion of |l|. For the 
stability analysis we transform the Cartesian coordinates 
of the particles to gyroscopic coordinates 



Xk + iVk = n ey^>(iVLt)[d k + f £], 
Xk-iyk = r&exp(- iQ£)[dk + £*.], 
z k = r b (z k + z* k ) 



(3) 



where Z k and £ fc are complex numbers defining the per- 
turbations to the circular orbit. In Eq. 10 the co- 
ordinates of the electron are defined by k = 1 and 
d\ ~ (7772 + m\)/m2 while the coordinates of the pro- 
ton are defined by k — 2 and d 2 ~ — (m2 + ml) /ml Q. 
Notice that the perturbations of Eq. ® are scaled by 
the light-cone distance tv The linear stability analysis 
proceeds like in Refs. 0, @, by substituting Eq. © 
into the equations of motion and expanding the equa- 
tions of motion to linear order. Analogously to Ref. Q, 
the tangent dynamics perpendicular to the plane (the z- 
direction) is decoupled from the tangent dynamics along 
the orbital plane to linear order. As found in Ref |7(, the 



linearized equations of motion have normal-mode solu- 
tions with every coordinate proportional to exp(Af2t/#). 
We shall see that |A| can be of the order of n or larger, 
such that Ail/ 9 is a frequency much larger than the or- 
bital frequency f2 (here the stiffness enters). Again a 
remarkable quasi-degeneracy is found between the pla- 
nar and the perpendicular normal- mode conditions |jj in 
the stiff-limit. If |A| is large enough (the stiff- limit), both 
conditions reduce to 



(^-)exp(-2A) = l, 



(4) 



where fj, is the reduced mass defined below Eq. Q and 
M = mi +7772 is the total mass of the atom. For hydrogen 
(fi/M) is a small factor of about (1/1824). Atomic orbits 
have 9 of the order of the fine structure constant such 
that the small parameter ^9 4 /M ~ 10~ 13 multiplying the 
exponential function in Eq.(0J determines that Re(A) = 

— a ~ — ln(v/ M//J,6 A ). For the first 13 excited states of 
hydrogen this a is in the interval 14.0 < \a\ < 18.0. The 
imaginary part of A can be an arbitrarily large multiple 
of 7r, such that the general solution to Eq. Q is 



A 



-a + nqi 



(5) 



where i = and q is an arbitrary integer. Notice that 
the real part of A is always negative, such that the tangent 
dynamics about circular orbits is stable in the stiff-limit 
with delay-only lj, at variance with the dynamics with 
advance and delay of Ref. [jj ■ 

The full normal-mode condition in the stiff-limit is de- 
rived in the same way of Ref. [jj and it will be published 
elsewhere. The perpendicular normal-mode condition is 
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and the planar normal-mode condition is 
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The coefficient of the term of order 1/A 2 is the main 
difference between Eqs. © and Q for 9 in the atomic 
range. The term on the right-hand side of Eqs. JOJ and 
(Q) with the 2/3 coefficient is due to the third derivative 
of the self- interaction force of Eq. J5J • The exact roots 
of Eqs. JJJ) and @ near the limiting-root © are defined 
respectively by 



K y {9) 
X z (9) 



-a xy + irqi + lei, 
-a z + nqi + ie 2 , 



(8) 
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where ei(9) and 62(6) are real and small numbers that 
depend on the orbit through 9. The stiff torus is formed 
from an initial circular orbit as follows; (i) the offending 
force against the velocity makes the electron loose radius 
by radiating energy on a slow timescale, and (ii) when the 
electron deviates enough from the circular orbit, the stiff 
nonlinear terms compensate the offending force of Eq. 
(J2J) and balance the real parts of Eq. (JHJ), the a's. The 
motion continues to be stiff because of the large imag- 
inary part of the A's of Eq. JHJ. In the following we 
outline a multiscale expansion about the balanced stiff 
torus to derive a resonance condition between the modes 
of tangent dynamics. The circular orbit recoils in a slow 
timescale to compensate the momentum of the radiated 
energy, such that the particle coordinates are described 
by a composition of a translation mode and a stiff mode, 
both with slowly varying amplitudes 



a very reasonable physical length. Notice also that the 
velocity associated with the stiff motion is given by 



v = —ilrba = 



irq 



(14) 



which is much lesser than the speed of light for 9 in the 
atomic scale. 

Assuming that |z*i| = \Ri\ = a, as estimated by Eq. 
(|12J) . we can use © to evaluate the angular momentum 
of the stiff torus along the orbital plane. Disregarding 
fast oscillatory terms, this angular momentum along the 
orbital plane is 

lx+Hy = Air2a 2 ^&iexp[t(e2-e? + 0)nt/0](15) 
9 

= ^Ja 2 exp[z(4-e? + 0)m/#], 



£ fc = A k (T)+u k (T)exp[(irqi + u$)Sn/6], 
Z k = B k (T) + R k (T)exp[(irqi + ieZ)nt/9}, 



(9) 



where the amplitudes u k and R k must be near the size a 
where the nonlinearity balances the negative real part of 
the linear modes. Because of the large frequency of the 
stiff motion, this balancing is established at a very small 
radius, as we estimate in the following. 

The correction to Eq. J2J is provided by the nonlinear 
terms of the Lorentz-Dirac force, as given in page 116 
of Ref 10] . Keeping only terms along the velocity, the 
self-interaction force along the electronic velocity is 



XljVi + 



(10) 



where we again multiplied by the convenient scaling fac- 
tor, as in Eq. (01. Using the stiff term of Eq. into 
Eq. to estimate the coefficient of the velocity term 
on the right-hand side of Eq. (|10fl yields 



(irqf a 2 



9} = 0, 



(11) 



where we used |ui| ~ |i?i| ~ a in Eq. ©, a being the 
radius of the stiff torus scaled by the light-cone radius 
n,. Equation Ijllfl estimates that along the stiff torus the 
stiff nonlinear term balances the offending force on the 
electron at the amplitude 
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ir 2 q 2 



(12) 



This distance is a small fraction of the orbital radius and 
it represent a length of 
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r b a 



9n 2 q 2 ' 



(13) 



which for 9 in the atomic scale is less than 30 classical 
electronic radia. This is the radius of the stiff torus, 



where we used formulas (0 for f2 and for on the sec- 
ond line of Eq. IjlSp. Notice that we chose the same 
q in Eq. 10 to cancel the fast oscillation in Eq. JTSJ. 
As illustrated in Fig. 1, the fast spinning angular mo- 
mentum should rotate with the circular orbit at the fre- 
quency ft, such that its phase should be equal to fit. This 
phase- velocity condition imposed on Eq. I|15|) yields the 
following resonance condition 







(16) 



The calculation of t\ and necessitates expanding the 
equations of motion to quadratic order to include the 
nonlinear stiff terms (the linear stability of the stiff 
torus). A purely harmonic solution to these nonlinear 
equations of tangent dynamics exists because the non- 
linear stiff terms introduce a correction proportional to 
a/9 into Eqs. 10 and 0, and the complete calcula- 
tion will be given elsewhere. We expect nevertheless the 
quasi-degeneracy to be preserved because the stiff limit 
depends basically on the time light takes to travel be- 
tween the particles. Once the stiff torus is very close to 
the circular orbit, t\ and should differ from ei and 62 
by a correction which is essentially the same because of 
the quasi-degeneracy property, plus a correction of order 
9, such that the difference can be well approximated by 



ei - £2 



b9 = b9 



(17) 



Condition (|16f) determines an orbital frequency propor- 
tional to a difference of two eigenvalues, a Rydberg-Ritz- 
like formula 



n = n(e 1 ~e 2 )/be = M -( ei - £2 ) 



(18) 



where we have used Kepler's law of Eq. 0. The exact 
calculation of b demands the inclusion of the nonlinear 
terms into Eqs. JfJJl and Q) and shall be given elsewhere. 



4 



Since b must be of order one, we henceforth set b = — 1 
into Eq. I|17(l as a qualitative approximation, yielding 

ei-e 2 + 6» = 0, (19) 

which has the solutions listed in Table 1 . Resonance l|16|) 
can also be derived from the solvability of the multiscale 
asymptotic solution with inclusion of the nonlinear terms 
by Freedholms's alternative^, and shall be given else- 
where 0| • In Ref ■ w e postulated heuristically a sim- 
ilar type of resonance for an analogous electromagnetic- 
like dynamical system, a resonance involving the real 
parts of the X's. The resonant orbits of Ref. were 
in the atomic scale and had the same qualitative QED- 
like dependence with the integer q, which seems to be a 
generic feature of these stiff electromagnetic resonances 
of the two-body problem. The root-searching problem of 
Eq. (|19|l is well posed and for each integer q it turns out 
that one can find a pair of the form (JSJ) if one sacrifices 9 
in Eqs. © and (0, i.e., 9 must be quantized! According 
to QED, circular Bohr orbits have maximal angular mo- 
menta and a radiative selection rule ( Al = ±1) restricts 
the decay from level k + 1 to level k only, i.e., circular 
orbits emit the first line of each spectroscopic series (Ly- 
man, Balmer, Ritz-Paschen, Brackett, etc.), the third 
column of Table 1. We have solved Eq. Ijl9(l and Eqs. 
© and (JZJ) with a Newton method in the complex A- 
plane. Every angular momentum 1/9 determined by Eq. 
(|19ll has the correct atomic magnitude for any q. The 
subset of resonant orbits listed in Table 1 has orbital fre- 
quencies near a circular line of QED with a few percent 
deviation. 
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Caption to Table 1: Numerically calculated angular 



momenta l z — 9~ x in units of e 2 /c, orbital frequencies 
in atomic units 137 3 2 (g 2 — ei), circular lines of QED in 
atomic units wqed = \ ~~ (fc+i) 5 ) ' an< ^ ^he values of 
the integer q of Eq. © . 

We identified the necessary existence of a fast timescale 
in the electromagnetic two-body motion near circular or- 
bits and explored the consequences of balancing the asso- 
ciated fast dynamics. The balancing of the fast dynam- 
ics predicts magnitudes in qualitative and quantitative 
agreement with QED. The stability analysis naturally in- 
volves an integer and the angular momenta are quantized. 
The emission lines are in reasonable quantitative agree- 
ment with those of QED and the angular momentum of 
the fast spinning motion is of the order of the lowest or- 
bital angular momentum, in agreement with QED where 
the electron has an angular momentum of spin given by 

Vsh/2 

Caption to Figure 1: Coulombian guiding-center cir- 
cular orbit of each particle, dashed lines. The trajectory 
of each particle is illustrated as the fast spiral encircling 
each guiding-center circular orbit (solid lines) . Arbitrary 
units, distances are not on scale, illustrative purposes 
only. 
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